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Abstract 

We address the removal of states with center-of-mass excitation from the SU(3) no- 
core shell model [SU(3)-NCSM] space, i.e., construction of the nonspurious subspace. 
A procedure is formulated based on solution of the null-space problem for the center- 
of-mass harmonic oscillator lowering operator -Bc°m\ operating at the level of SU(3) 
irreducible representations. Isolation of the center-of-mass free subspace for the SU(3)- 
NCSM provides the foundation for exact removal of center-of-mass dynamics in the 
proposed Sp(3, R) symplectic no-core shell model. We outline the construction process for 
the matrix representation of Bc°m. , present the algorithm for obtaining the nonspurious 
space, and examine the dimensions obtained for center-of-mass free SU(3) subspaces in 
representative light nuclei. 
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1. Introduction 

The ability to carry out no-core configuration interaction calculations of light nuclei, 
in the no-core shell model (NCSM) [1-3], has made a significant contribution to recent 
progress in the ab initio description of nuclei. However, the dimensionality of the nuclear 
model space becomes computationally prohibitive as the number of active nucleons and 
orbitals increases. Symmetry can play a significant role in addressing this problem, by 
assisting in the selection of the physically relevant portions of the model space. Elliott [4- 
6] explored the U(3) D SU(3) symmetry of the harmonic oscillator, which serves as 
an organizational scheme for quadrupole deformation and rotation in the nuclear shell 
model. The U(3) algebra of the oscillator is contained in a larger Sp(3,R) algebra, which 
is found to have a close connection both to the dominant components of the nuclear 
Hamiltonian and to nuclear collective motion [7-9] . Following the discovery of evidence 
for SU(3) and Sp(3,M) symmetry in the nuclear eigenstates obtained in conventional 
NCSM calculations for light nuclei [10, 11], the NCSM has been reformulated in terms 
of an SU(3)-based model space, in the SU(3) no-core shell model [SU(3)-NCSM] [12, 13]. 
This development provides the foundations for future realization of a symplectic no-core 
shell model (Sp-NCSM) [14], making full use of Sp(3,R) symmetry. 

In predicting physical properties of the nucleus, only the intrinsic dynamics of the 
nucleus is of interest, not the center-of-mass motion. In principle, the center-of-mass 
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motion may be eliminated from the problem by explicitly changing variables to relative 
coordinates. However, the nuclear many-body state must be antisymmctrized, and with 
increasing nucleon number this process rapidly becomes prohibitive in relative coordi- 
nates [3]. On the other hand, for a many-body basis constructed from antisymmctrized 
products of single-particle states, as in the NCSM, antisymmetrization is straightfor- 
ward, but both center-of-mass (spurious) and intrinsic excitations are included in the 
model space [15, 16]. Consequently many of the nuclear eigenstates obtained by diag- 
onalizing the Hamiltonian will carry center-of-mass excitation. Moreover, exact factor- 
ization of the many-body wave function into center-of-mass and intrinsic parts — i.e., 
ip{ri](Ji) = ipc.m.{R)' l Pin(rij;<Ti), where V'c.m.(-R) depends on the center-of-mass coordi- 
nate, and ipin{fij;iJi) depends on relative coordinates and intrinsic spin degrees of 
freedom — is obtained only within certain specific truncations of the model space. In 
the context of the conventional NCSM, for which the basis is constructed from antisym- 
mctrized products of oscillator wave functions, factorization is obtained if truncation is 
carried out according to the iV max scheme, that is, by the total number of oscillator quanta 
for the many-body state [3]. Within the SU(3)-NCSM, factorization is also retained in 
model spaces which have furthermore been truncated according to SU(3) and spin quan- 
tum numbers (see Section 3). Otherwise, factorization of intrinsic and center-of-mass 
wave functions is in general incomplete, in which case all eigenstates contain some spu- 
rious contribution, presenting challenges for the study of the intrinsic excitations {e.g., 
Refs. [17-19]). 

The usual approach to addressing the problem of spurious states consists of modifying 
the nuclear Hamiltonian, working within a model space which supports factorization, by 
adding a Lawson term [20] proportional the harmonic oscillator Hamiltonian i/ c . m . for 
the center-of-mass coordinate. After diagonalization, the spurious states remain in the 
spectrum but are shifted to an excitation energy above the low-lying intrinsic states of 
interest. However, the other possibility is to modify the model space, so as to explicitly 
remove the entire space of spurious states, before the Hamiltonian is diagonalized. This 
approach may be made feasible through the use of SU(3) symmetry [21-25]. After the 
spurious states are removed, the remaining center-of-mass free (CMF) states form a 
suitable model space for describing the intrinsic dynamics of the nucleus. 

The CMF states are associated with the lowest eigenvalue of the center-of-mass Hamil- 
tonian i/ c . m . and thus possess a harmonic oscillator 0s wave function in the center-of-mass 
degrees of freedom. They also therefore constitute the null space of the center-of-mass 
harmonic oscillator lowering operator £?c°m.- The approach developed here for isolating 
the CMF space is based on solution of the null-space problem for this operator, formu- 
lated at the level of subspaces of definite SU(3) symmetry and intrinsic spin. 

A principal motivation is to enable exact separation of center-of-mass dynamics in 
extensions to the SU(3)-NCSM, in particular the Sp-NCSM. The many-body basis for 
the SU(3)-NCSM is the starting point for defining the basis for the Sp-NCSM [14]. 
Briefly, an Sp(3,R) irreducible representation (irrcp) is constructed from an extremal 
SU(3) state [14], i.e., possessing the fewest oscillator quanta. Other states within the 
irrep are constructed by acting on this extremal state with an SU(3)-coupled product 
[^(20) x 4(20) _ x j4 (20)](A M ) ) w hcrc A^ 20 ) is the translationally-invariant Sp(3,R) raising 
operator. Each application of A^ 20 ' contributes two oscillator quanta to the system. 
Since the action of A^ 20 ) does not introduce center-of-mass excitation, an Sp(3,R) irrep 
is free of center-of-mass excitation as long as it is built from an extremal SU(3) state 
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which is free of such excitation. The present results therefore provide the foundation for 
obtaining the CMF model space for the Sp-NCSM. 

The construction of the SU(3)-NCSM basis, before removal of spurious contributions, 
is first outlined (Section 2). We then introduce an algorithm for identification of the 
CMF subspace within the SU(3)-NCSM model space, based on constructing the matrix 
representation of i?c°m. between SU(3)-coupled subspaces and solving for its null space 
(Section 3). Finally, we summarize the dimensions obtained for CMF SU(3) subspaces 
in representative light nuclei (Section 4). Preliminary results were reported in Rcf. [26]. 

2. SU(3)-NCSM basis states 

The many-body basis states for the SU(3)-NCSM have good SU(3) x SU(2) quantum 
numbers, where the SU(3) symmetry label (A/i) characterizes the spatial degrees of free- 
dom, according to the Elliott classification, while the SU(2) label S describes the total 
intrinsic spin angular momentum. The creation operator for a nucleon in a given major 
oscillator shell rj comprises an SU(3) x SU(2) tensor oLp)i/2i where the labels denote 
(A/i)5 = (770)5 [27]. The operators 0^0)1/2 ar e then used as the fundamental units in 
building up an SU(3) x SU(2)-coupled nuclear state. 

Specifically, each SU(3)-NCSM basis state is characterized by a definite distribution 
of nucleons over the major shells. First, all nucleons in the each major shell 77 are 
combined to form a configuration of the type [oLq^^ x a ( I( o)i/2 x ' ' •]( A ' ,A^7, )' S ' , , with 
good SU(3) x SU(2) coupling [28], separately for the protons and neutrons. Then, such 
configurations for individual major shells are coupled successively to form a total proton 
state and total neutron state carrying good SU(3) x SU(2) quantum numbers. Finally, the 
proton state and neutron state are coupled, to give the SU(3)-NCSM basis state. Since 
the major shells have definite occupations, the state may be classified, as usual in the 
oscillator basis for the NCSM [3], by the number iV cx of harmonic oscillator excitations, 
taken relative to the minimal number of oscillator quanta possible for the given number 
of protons and neutrons. 

The resulting SU(3)-NCSM basis state has the form, with all coupling labels shown 
explicitly, 

(l((((7 P ,o x 7 P a) Ppi0 " P '° x Ipfl) 9 " 1 "" 1 x 7p,3) Pp ' 2 ^ 2 - x 7 P ,r, mil J pp " p > 

X |(((( 7 n,0 X 7n,l) P "'°^° X ln ^^ X ln ^^ ... X 7 n,„ ma J P " W ">) (1) 

Here the symbol -y Stn represents the labels 7 = [/1, /2, . . . , /jv]a(A/x)S needed to com- 
pletely specify the coupling of nucleons of type s (i.e., protons or neutrons) within ma- 
jor shell 77 [28]. Specifically, each major shell has associated with it a U(iV) algebra 
[N = (77 + 1) (77 + 2)] consisting of bilinears of creation and annihilation operators, for 
which the irreps are labeled by [/1, /2, ■ • ■ , /jv], where we consider only antisymmetric ir- 
reps, and J2i=i /» equals the occupation of the shell. Within a U(iV) irrep, a multiplicity 
index a is required to distinguish SU(3) x SU(2) irreps with the same quantum numbers 
(Xfi)S, yielding the labeling scheme 

U(iV) D SU(3) x SU(2) . (2) 

[/1,/2,-,/iv] a (Ap) S 
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The symbols uj s , v in (1) then indicate the SU(3) x SU(2) coupling labels [uj = (Xp)S] of 
successive shells, and p StTj denotes the multiplicity index for this coupling. Finally we 
have total couplings uj p = (X p p p )S p for the protons, uj n = (X n p p )S n for the neutrons, 
and co = (Xp)S for the entire basis state, with corresponding multiplicity indices p p , p n , 
and p, respectively. 

The expression in (1) represents not just a single state but an entire set of states, 
with various values for the quantum numbers associated with the branching of SU(3) x 
SU(2) into angular momentum subalgebras: the orbital (spatial) angular momentum L, 
the inner multiplicity k for this L within the SU(3) irrcp (Xp), and the total angular 
momentum J, as well as its z-projection M. However, these states share the same 
"internal" microscopic structure, given by the same couplings of the particles at the level 
of SU(3) x SU(2). Therefore, they may be thought of as a single reduced state for certain 
purposes, in particular, evaluation of reduced matrix elements under the SU(3) x SU(2) 
Wigner-Eckart theorem and, as we shall see, identification of CMF linear combinations. 
The analogy is to angular momentum theory, where one may consider the states \JM), 
for different M, to be substates of a single state || J), more formally, a tensorial set or 
SU(2) irrep. 

3. Construction of the CMF subspace 

The separation of the many-body space into CMF and spurious parts simplifies in the 
context of an SU(3) x SU(2)-coupled basis, since the process may be carried out inde- 
pendently within subspaces characterized by definite SU(3) x SU(2) quantum numbers. 
To start with, the center-of-mass Hamiltonian i? c . m . does not connect states involving 
different numbers iV cx of oscillator excitations, i.e., [H c , m , 7 N c ^\ = 0. It is due to this 
property that the usual N max truncation scheme for the NCSM [3] permits an exact fac- 
torization of center-of-mass and intrinsic wave functions. 1 This property also implies that 
the separation of CMF states may be carried out separately within the space of states 
with each specific number of oscillator quanta, which we denote by Wjv ex • Furthermore, 
-Hc.m. commutes with the SU(3) generators and may therefore be diagonalized within a 
subspace with good SU(3) quantum numbers [22]. As an operator acting only on spatial 
degrees of freedom, ff c . m . also commutes with the total spin operators for protons and 
neutrons, as well as their combined spin operator. Thus, collecting these properties, the 
separation of CMF states may be carried out separately within subspaces of given iV ox , 
(Xp), S p , S n , and 5, which we denote by WN ox [(Xp)S p S n S]. 

In considering how to extract CMF states, we note that these states are defined by 
their relation to the center-of-mass harmonic oscillator raising and lowering operators. 
The center-of-mass raising operator Bl. m ., which is an L = 1 operator, furthermore 
constitutes an SU(3) (10) tensor, with components m° l=i m . This operator may be 
written in terms of single-particle harmonic-oscillator raising operators b^ 10 ^ as [14, 27] 

« = 4fl> t(10) «- ( 3 ) 



lr The 7V max truncation is to many-body states with < N cx < 7V max . For NCSM calculations with 
parity-conserving interactions, the two parity subspaces, obtained for 7V ex even or odd, respectively, may 
furthermore be considered separately. 
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The corresponding lowering operator -Bc.m. similarly has the form 

B^ = ^=J2~b^(i). (4) 

The center-of-mass harmonic oscillator Hamiltonian is built from these operators as 

H c . m . = ^(Bt^.B^ + l), (5) 

where the dot indicates a spherical tensor scalar product. The CMF states are then 
defined by the property that they have no center-of-mass excitations, i.e., they have zero 
eigenvalue for the center-of-mass number operator N c , m , = ■ £?c°m.- Equivalently, 

however, they arc identified by the property that they are annihilated by the center-of- 
mass lowering operator Bc°m.- 

Either of these criteria allow the problem of identifying CMF states to be for- 
mulated as a null-space problem. We seek the subspace W^- MF [(Xfi)S p S n S] of 
W Nex [(X/j^SpSnS] consisting of states |* CMF ) such that iV c . m .|* CMF ) = or, equiva- 
lently, B^ c ^ L=1 M \^ CMF ) = 0. In practice, this means first representing the operator as 
a matrix with respect to the SU(3)-NCSM basis and then solving for a complete set of 
null vectors of this matrix. The result yields new basis states, for W^ MF [(X/j,)S p S n S], as 
linear combinations of the original basis states, for Wjv ex [(Xjj,)S p S n S]. 

Although we could in principle search for the null space of either N c , m , or 5c°m., there 
is a practical advantage to working with _Bc°m.- While iV c . m . is a two-body operator, £?c°m. 
is simply a one-body operator. Consequently, evaluation of matrix elements is computa- 
tionally less involved. Note that iV c . m . acts within the space Wjy ex , i.e., conserving the 
number of oscillator excitations, but Bc°m. connects the space Wjv 6x to the next lower 
space Wjv ex -i, which is significantly smaller in dimension than Wjv«- One might there- 
fore expect the null space problem for _Bc°m. to be of lower dimensionality than that for 
N c , m ,, namely, involving a matrix of the same column dimension (~ dim W0v ex ) but lower 
row dimension (<~ dim Wjv ex - 1 ) • However, this simple relation of dimensions is modified 
once SU(3) selection rules are considered (see below), and in practice the difference in 
dimensionality of the two problems is largely eliminated. 

The problem of identifying CMF states is simplified by the realization that it may be 
formulated entirely at the level of reduced states. Recall that i/ c . m . commutes with the 
SU(3) generators, which connect states within an irrep. Consequently, an SU(3) x SU(2)- 
reduced state is CMF if and only if its substates, labeled by nLJ (and M), are all CMF. 
We need thus only find a basis of SU(3) x SU "(2)- reduced states, for W^ F [(Xfi)S p S n S], 
independent of kLJM. We also need only consider the SU(3) x SU {2)- reduced matrix 
elements of the operator Sc°m., rather than the matrix elements among individual kLJM 
substates, which are much greater in number. 

To explicitly relate the null space problem at the level of individual states to that for 
SU(3) x SU(2)-reduced matrix elements, let us return for a moment to kLJM states and 
observe that any CMF state \^%^J[{X^)S p S n S; kLJM]) within the W Nex [(Xfj,)S p S n S] 
subspace must satisfy 

^ Nm -i[(X'^)S p S' n S'; K 'L'J'M']\Bi ^ = 0, (6) 
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for every state \^ Nex -i[(X' fi')S p S' n S'; k'L'J'M']) e Wn^-i- Since we are working with 
states of good angular momentum, we can immediately rewrite the condition in terms of 
a reduced matrix element as 

(* Wex _ 1 [(AV0^^5^'i'^ll^m.,L =ill*^ F [(AM)5 P 5 ri 5;KLJ]) = 0. (7) 

Note that, since _Bc°m. acts only on spatial degrees of freedom, we actually need only 
consider the case (S p S' n S') — (S p S n S). The ordinary SU(2)-reduced matrix element 
in (7) is related to the SU(3) x SU(2)-rcduced matrix element of B ( c °^. by the SU(3) 
Wigner-Eckart theorem (and LS'-coupling relations) , as 



_ ( .j+L'+i+sf.pU j s\ /(am) (oi) (Ay; 

_ [ > L J \ J' V 1 )\ kL 1 k'V 

x (^.^(AyjSpS^llB^II^KA/iJS^nS]), (8) 



where J = (2J + l) 1 / 2 , and the quantity in parentheses is an SU(3) Clebsch-Gordan 
coefficient [29]. The condition that the complete set of SU(2)-reduced matrix elements 
appearing on the left-hand side of (8) vanish is equivalent to the condition that the single 
SU(3) x SU(2)-reduced matrix element on the right-hand side vanish. 

Thus, for the CMF states, we seek SU(3) x SU(2)-reduced states ||*™ f [(Am); S p S n S\) 
such that 

(^-iKAVOSpSnSlllB^II^MSpSnS]) - 0, (9) 

for all possible reduced states ||* JVex -i[(Ay) < S' p S' rl S']) e WV ex -i [(X' n')S p S n S]. The sub- 
spaces Wjv ex -i[(A' jj,')S p S n S] which may be linked with WN BX [(Xfj,)S p S n S] through £?c°m. 
are restricted by the SU(3) tensor character of £?c°m.- Specifically, (X'jj, 1 ) must be con- 
tained in the product (A//) x (01) which, from the general rules of SU(3) coupling [30], 
may be seen to consist of 



(Am) x (oi) = < 



'(01) A = M = 

(0/i+l)®(l/i-l) A = 0,a*>1 

(Al)e(A-lO) A>1,m = 

k (A M + i)e (A + l m- i)e (A- 1 m) A>1,m>1- 



(10) 



The problem of finding reduced states which satisfy (9) can be converted into 
searching for the null space of a matrix, the entries of which are the SU(3) x SU(2)- 
reduced matrix elements of -Bc°m. between the SU(3) x SU(2)-reduced basis states for the 
WN ox [(Xn)S p S n S] subspace and the SU(3) x SU(2)-reduced basis states for each of the 
possible subspaces Wjv ex -i[(A' fi')S p S n S]. The resulting matrix has the form illustrated 
in Fig. 1, where the horizontal dashed lines delimit submatrices corresponding to the dif- 
ferent final spaces with (A'//) = (A^m'i), (A2M2); • • ■> as allowed by the selection rule (10). 
(There will be at most three such submatrices.) The entries of the null vectors then give 
the expansion coefficients for the basis states for the CMF space W™ F [(X[i)S p S n S] in 
terms of the original basis states of WN ex [(Xfj,)S p S n S]. 

The SU(3)-reduced matrix elements of £?c°m. entering into the matrix of Fig. 1 can be 
calculated numerically using methods from Ref. [31]. In particular, the SU(3) x SU(2) 
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^iVex-l.lll-^C.m.H^JVex.l/ 

||s (01) ll* (v) ) 

^ATex-l^ll-^c.m.H^jv^,!/ 



\*AT ex -l,lll- D c.m.||*^ ex; i/ 

/, T/ ( A 2A4) II 5(01) ||,T,(V) \ 



ii^(Ol) v 
/.TvW^i) || r( 01 )||iT;( A ^) \ 



(^l.lll^ll^i) 



Figure 1: Form of the matrix representation of 6j m , as reduced matrix elements between the SU(3)- 
NCSM basis states for the Wjv ex [(Xfi)S p S n S] subspace (identified with the columns) and the SU(3)- 
NCSM basis states for each of the possible subspaces VKjv cx -i [(X' fi')S p S n S] (identified with the rows). 



matrix elements of a one-body operator such as Bc°m. (or an n-body operator in general) 
can readily be computed, once the operator is expressed in second-quantized form in 



t 

(770)1/2, 



terms of SU(3) x SU(2)-coupled products of creation operators a 

operators fi( 0n )i/2 [27]. For the center-of-mass annihilation operator\Bc U m., the second- 
quantized form is obtained as 



and annihilation 

(01) 



B(oi) 



1 

71 



E 



(77 + l)(r7 + 2) 



(V\\b m \\v + 1) 



I (r,0)l/2 X a (0,r,+ l)l/2 



(01)0 



(11) 



where (?7||&( 01 )||?7 + 1) = 



4. Dimensions of CMF spaces 

Let us now examine the dimensions of the CMF subspaces obtained by the methods 
of Section 3. Our primary interest is in the distribution of CMF and spurious states 
with respect to the SU(3) x SU(2) quantum numbers. The Sp-NCSM approach requires 
identification of these CMF SU(3) x SU(2) irreps, at low N cx , from which Sp(3,R) irreps 
extending to high N cx are created by repeated action of the symplectic raising operator 
(Section 1). 

To begin with, the simplest illustration we might consider is the model space for the 
deuteron ( 2 H), which is shown along with that for the triton ( 3 H) in Fig. 2. The full 
SU(3)-NCSM space for the deuteron may be broken into subspaces WN ex [(\/j,)^S], with 
5 = and 1. The dimensions of these subspaces — by which we specifically mean the 
number of SU(3) x SU(2)-reduced basis states, not the total number of kLJ basis states, 
which would be much higher — are indicated by the areas of the outer, light circles in 
Fig. 2. Then, for each value of N ex , it is found that solution of the null space problem 
for the center-of-mass annihilation operator yields two CMF reduced states. Both have 
(Ait) = (iV ex 0), one with S = and one with S = 1. These CMF subspaces are indicated 
by the inner, dark circles in Fig. 2. 
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Figure 2: Dimensions of the SU(3) X SU(2)-reduced subspaces of the SU(3)-NCSM model space for 
(a) the deuteron ( 2 H) and (b) the triton ( 3 H). The area of the light (outer) circle indicates the full 
dimension, i.e., of Wpf ox [(Xii)SpS n S], and the area of the dark (inner) circle indicates the dimension of 
its CMF subspace, i.e., of W^ F [(\fj.)S p S n S]. Subspaces are arranged by 7V CX (labeled by N cx hQ at 
top), SU(3) quantum numbers (at bottom), and spin quantum numbers (S p S n S) = (^iS) (at left). 



The deuteron provides a particularly illuminating example, since the quantum num- 
bers of the CMF spaces obtained in Fig. 2 may be understood through simple arguments. 
Let iVc.m. and N le \ denote the number of oscillator excitations of the center-of-mass and 
relative degrees of freedom, respectively, so iV cx = -/V c m . + N rs \. The CMF condition 
imposes iV c . m . = and thus N rc \ = N cx . For the two-particle system, the transforma- 
tion between single-particle and relative coordinates is straightforward. There is only a 
single relative coordinate vector, and the harmonic oscillator in this coordinate is equiv- 
alent to the harmonic oscillator for a single particle in three dimensions. For a given 
-/V cx , it will therefore carry SU(3) quantum numbers (7V ex 0). The center-of-mass os- 
cillator carries (00) for a CMF state. Therefore the CMF state as a whole will have 
(Xfi) — (N CK 0) x (00) = (7V OX 0), as well. Since the deuteron consists of distinguishable 
particles, the coupling of spins is independent of the spatial degrees of freedom, and both 
5 = and 1 are obtained. 

In general, as we move beyond the two-body problem, the fraction of the model space 
dimension which corresponds to CMF states, versus that which corresponds to spurious 
states, depends both on the number of nucleons and on the number iV cx of oscillator 
quanta. The spurious contribution at N cx = is identically vanishing, as shown by 
Elliott and Skyrme [15]. For a given nucleus, the spurious contribution becomes an 
increasing fraction of the total space with increasing iV ex , while, for a given level of 
excitation N ex , the spurious contribution becomes a less significant fraction of the total 
space with increasing atomic mass A. These trends are already in evidence in comparing 
deuterium with tritium (Fig. 2) but may be seen more systematically in Fig. 3(a), where 
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Figure 3: Dimension of the SU(3)-NCSM model space, decomposed according to N ex : (a) by number 
of SU(3) X SU(2)-reduced basis states, (b) after branching to J = states, and (c) after branching to 
J = 2 states. Both full (solid curves) and CMF (dashed curves) dimensions are shown, for 2 H, 4 He, 6 Li, 
8 Be, and 12 C. 

the dimensions of the full spaces (solid curves) and CMF spaces (dashed curves) are 
shown, as functions of iV cx , for several nuclei with 2 < A < 12. 

Qualitatively, the pattern in the evolution of the CMF fraction may be understood 
from simple counting arguments, by considering the possible ways of distributing iV cx 
oscillator quanta over the 3^4 oscillator degrees of freedom of the nuclear system. These 
may be decomposed into three center-of-mass oscillator degrees of freedom and 3^4 — 3 
relative oscillator degrees of freedom. If many oscillator quanta are to be distributed over 
few degrees of freedom (N cx 3> A), most of the distributions will allocate at least one 
oscillator quantum to the center-of-mass degrees of freedom, leading to a high proportion 
of spurious states, i.e., a smaller CMF fraction. Conversely, if few oscillator quanta are to 
be distributed over many degrees of freedom (iV cx <C A), most distributions will "miss" 
these three center-of-mass degrees of freedom, i.e., not allocate any oscillator quanta to 
them, leading to a low proportion of spurious states, i.e., a larger CMF fraction. 2 

We focus now on the detailed distribution of CMF and spurious states over the SU(3) x 
SU(2) subspaces of Wjv ex , that is, the N e x(\/J,)S p S n S subspaces of Section 3. First, it 
should be noted that the dimensions of these subspaces vary widely. The dimensions of 
the SU(3) xSU(2) subspaces for |He 2 , |Li 3 , and |Be 4 , with iV cx < 8, 6, and 4, respectively, 
are shown in Fig. 4 — these dimensions vary from one to 640. To accommodate this 
range, the circles in Fig. 4 have been sized according to a power law scale, such that 
a doubling in radius represents a tenfold increase in dimension. However, any such 
nonlinear rescaling precludes meaningful visual comparison of the dimensions of the full 
and CMF subspaces, as was possible in Fig. 2. Consequently, we examine the CMF 



2 Quantitatively, such counting arguments give an estimated CMF fraction (dim PV^ MF )/(dim Wjv ex ) 
Ri (3A — 1)-/(7V CX + 3A — 1)-, where m— = m(m — 1) • ■ ■ (m — n + 1). This expression gives the exact 
result for the dcuteron, when applied to the total space, i.e., counting all M substates, but is of only 
approximate validity for other nuclei, due to neglect of antisymmetrization. For large 7V CX , we obtain a 
CMF fraction falling as N~ x 3 . 

9 



OfiQ. 2MI 3m Am 



5m 



6m 



im 



%m 



(1,1,2) 
Co (1,1,1) 

^ (1,0,1) 
(0,1,1) 
W (0,0,0) 



(a);; 
4 He[ 



• • • . • 



• • • • • 



• •••••i 

•••••• 



••••••• 

••••••• 

••••••• 

••••••• 



© ^ CN ^ <0 to CN ^ <0 <0 ~H en <0 CN ^ <0 <0 CN ^ CO <0 CN ^ O r ~H <o CN ^ CO <0 CN ^ O r 

o~ o~ cn~ o~ rn o~ ri o~ ri ^ m" <rT o~ o~ c-T ri *o ^-T o~ cn co »n r-" o~ cn ~h co ©~ cn >n oo~ 



ofin imi 2*n 



3W1 



4^n 



sm 



6m 



(3/2,3/2,3) 
(3/2,3/2,2) 
(3/2,3/2,1) 
££(3/2,3/2,0) 
=(3/2,1/2,2) 
"=,(3/2,1/2,1) 
Co(l/2,3/2,2) 
(1/2,3/2,1) 
(1/2,1/2,1) 
(1/2,1/2,0) 



(2,2,4) 
(2,2,3) 
(2,2,2) 
(2,2,1) 
(2,2,0) 
(2,1,3) 
(2,1,2) 
(2,1,1) 
(1,2,3) 
(1,2,2) 



(b) 
6 Li 



• • • • 



• • • • • 



o O O O O o 



•••••• 



• •••••••• • 

• •••••••• • 

>•••••••• • 

>••••••••• • 

!•••••••••• 

• •••••• 

• •••••• ••••••••••• 

• •••••• ••••••••••• 

•••••••••••••••••• 

•••••••••••••••••a 



— o 

o"<N 



°„"~'., .. ., 



CN O 
o"cN 



O — i CN ^ 
O — " O cn CN ^ 



„ °i ^"J, ^ „ °i „ ^"J, . °i ^"i ^ ^ „ ~t. ~ 

~h o~ CN~ — r cn <rT I"-" CN~ ~h cn o~ CN >n 00 



om 



im 



2m 



3m 



4m 



CO 

c 
CO 



CO 



(1,2,1) 



(2,0,2) 
(1,1,2) 
(1,1,1) 
(1,1,0) 
(0,2,2) 
(1,0,1) 
(0,1,1) 
(0,0,0) 



(c) 
8 Be 



°„ 

— " o" CN 































o 


• • 










































o 


• • 










































o 


o 








































■ 


o 


• • 




















































































o 


o 


• • 


o 


o 


o 


o 
















o 


o 


o 




o 








° 


o 


• • 


o 


o 


o 



















o 




o 












o 


o 


• • 




o 


o 


o 
































• 


o 


• • 


o 


o 


o 


o 
















o 




o 




o 








o 


o 


• • 


o 


o 


o 


o 
















o 




o 












o 


o 


• • 


o 


o 


o 


o 
































o 


o 


• • 


o 


o 












I • • 







o 







o 


o 


G 


o 





o 




O 


o 


oo 


O 


• • 


o 


O 





1 • • 


o 


o 


o 


• 


o 


o 


o 


o 


o 


o 


o 




o 


o 


oo 


o 


o 


o 


o 


o 


o 


1 • • 




o 


o 


• 


o 


o 


O 


o 


o 


o 


o 




o 


o 


• • 


o 


• • 


o 


o 


o 














o 


o 




o 








o 


o 


• • 







o 


o 


o 




1 • • 


O 


o 


o 


• 


O 


O 





o 


o 


o 


o 




O 


o 


• • 


o 


o 


O 


o 


o 


o 


1 • • 


o 





o 


• 


O 


O 


o 


o 


o 


o 







O 


o 


• • 


o 


o 


O 


o 





o 


1 • * 

J 1 l_ 


o 

l_ 


, l_ 


o 

1 


• • 

1 l_ 


O 

i_ 




1_ 


o 

1 


o 

1 




1 


o 

1 


o 


o 

l_ 


o 

l_ 


o 

l_ 


• • 

1 1 


o 

_l_ 


o 

1_ 


O 

i_ 


o 

1_ 


o 

1 




i 


°„ ^ 






CI 


— ^ <o 




<N 




co 


o 


CN 




c 




q 


<N — ^ 




cn 


o_ 


EN 







Figure 4: Dimensions of the SU(3) X SU(2)-reduced subspaces for (a) 4 He, (b) 6 Li, and (c) 8 Be. The 
size of each circle represents the dimension, scaled such that a doubling in radius represents a tenfold 
increase in dimension. Subspaces are arranged by A r ex (A/^)5 p 5 n S as indicated in the caption to Fig. 2. 
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Figure 5: Ratio of CMF dimension to total dimension, for the SU(3) X SU(2)-reduced subspaces for 
(a) 4 He, (b) 6 Li, and (c) 8 Be. The area of the darkened portion of each square, relative to the total 
area, indicates the dimension of the CMF subspace, i.e., of W^- MF [(\fj.)SpS n S] relative to that of the 
full subspace, i.e., of Wjv ox [(XlJ,)SpS n S]. Subspaces are arranged by N cx (\fi)S p S n S as indicated in the 
caption to Fig. 2. 
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fractions, i.e., the ratio (di-mW™ F [(\fi)SpS n S])/(dimWN cx [{X^)SpS n S}) 7 separately in 
Fig. 5, for these same nuclei. 

As already noted, the CMF fractions decrease with increasing N cx and approach unity, 
for a given iV cx , with increasing mass. However, within each 7V 0X space, demarcated by 
vertical dashed lines in Fig. 5, a further dominant trend may be noticed in the variation of 
CMF fractions among the SU(3) x SU(2) subspaces. It is observed that the spurious states 
are preferentially found in the subspaces with SU(3) quantum numbers corresponding 
to the lowest values of the SU(3) second-order Casimir invariant C2, while the CMF 
states are preferentially found in the subspaces corresponding to the highest eigenvalues 
of C2. In Fig. 5 (as well as in Figs. 2 and 4), within each VFjv ox , the SU(3) irrep labels 
(A/i) are ordered from left to right by increasing eigenvalue of C'2 (given by (C2) = 
A 2 + A/i + /j 2 + 3A + 3/i), where labels which are degenerate with respect to C2 are then 
ordered by increasing A. The lowest CMF fractions are (predominantly) found at left, 
and the highest (predominantly) at right. While some evidence may be seen for patterns 
in the distribution with respect to the spin labels, these are not as clear or consistent. 

So far, we have considered dimensions at the level of SU(3) x SU(2)-reduced basis 
states. We may also deduce from these the full and CMF dimensions in terms of J-couplcd 
states. Although much of the computational process in an SU(3)-NCSM calculation 
can be carried out in terms of SU(3) x SU(2)-reduced states, and it is these reduced 
states which are relevant to the definition of the Sp-NCSM model space, the Hamiltonian 
matrix in an SU(3)-NCSM calculation must ultimately be realized in terms of basis 
states of definite angular momentum J. As outlined in Section 2, the J states are 
obtained by first branching each SU(3) x SU(2)-reduced state (1) to states of definite 
orbital angular momentum L. The L values contained within an SU(3) irrep (A/z) are 
given by the SU(3) — > SO(3) branching rule [4]. 3 Then L is coupled with S to yield 
states of total angular momentum J according to the usual coupling rules for angular 
momentum addition. The dimensions of the resulting J-spaces — both the full space 
and its CMF portion — are shown for several light nuclei in Fig. 3(b,c), for J = and 
J = 2, respectively. Although the dimensions of the J-spaces are calculated here via 
the SU(3) x SU(2) coupling scheme, the results obtained are generally applicable to any 
J-coupled scheme for the NCSM, in an N max truncation. 

5. Conclusions 

The separation or elimination of spurious center-of-mass excitations is essential to 
the problem of determining the intrinsic structure of nuclei. If an SU(3)-coupled har- 
monic oscillator basis is used for the many-body problem, as in the SU(3)-NCSM, the 
separation may be carried out at the level of SU(3) irreps, in particular, within sub- 
spaces of fixed number of oscillator quanta and (X/j,)S p S n S labels. We have formulated 
the problem of finding the CMF subspace as a matrix null-space problem, based on 
the SU(3) x SU(2)-reduced matrix elements of the center-of-mass annihilation operator, 
which is solved independently for each SU(3) x SU(2) subspace of the full model space. It 



3 For each value K = min(A, fi) mod 2, . . . , min(A, n) — 2, min(A, fi), states are obtained with L = 
K, K + 1, . . . , _fsT + max(A, fi), with the exception that L = max(A, fi) mod 2, . . . , max(A, /1) — 2, max(A, fi) 
if K = 0. 
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is therefore possible to remove spurious contributions from the SU(3)-NCSM model space 
prior to diagonalization of the Hamiltonian, rather than through a Lawson term. In the 
context of the SU(3)-NCSM per se, this raises the possibility of substantial reductions in 
dimensionality of the problem, principally in the high N cx subspaces of lighter nuclei. 

However, a more essential application lies in providing the foundation for ensuring 
exact separation of center-of-mass and intrinsic dynamics, or removal of spurious con- 
tributions, in the Sp-NCSM. The purpose of the Sp-NCSM is to incorporate physically 
relevant portions of the model space extending to much higher numbers of oscillator 
quanta, beyond those which can be practically reached if one must retain the complete 
N ex (Xii)S p S n S subspaces of the SU(3)-NCSM. However, factorization of ccntcr-of-mass 
and intrinsic wave functions is still guaranteed if the Sp(3, R) basis is constructed starting 
from SU(3)-NCSM extremal states which are free of spurious excitation. The present ap- 
proach may also serve as the starting point for eliminating spurious admixtures in other 
extensions to the SU(3)-NCSM, such as an adaptation of the importance truncation 
scheme [32] to the SU(3)-NCSM. 
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